Hydrodynamics of topological defects in nematic liquid crystals 
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We show that back-flow, the coupling between the order parameter and the velocity fields, has a 
significant effect on the motion of defects in nematic liquid crystals. In particular the defect speed 
can depend strongly on the topological strength in two dimensions and on the sense of rotation of 
the director about the core in three dimensions. 
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Topological defects arise in all areas of physics from 
cosmic strings to vortices in superfluid helium . Al- 
though the physical systems are very different, many as- 
pects of the observed phenomena match even quantita- 
tively, making it possible to test cosmological predictions 
experimentally in condensed matter systems De- 
fects are classified according to their topological strength, 
and in many cases, the symmetries of the field equations 
lead to dynamics where defects of opposite topological 
strength can be mapped into each other. In this Letter 
we show an example, of topological line defects in liquid 
crystals where this symmetry is broken, due to the 
coupling to an additional field, the flow-field, or to elastic 
constants that arc not equal. 

In liquid crystals the topological defects are moving 
within a liquid and therefore one must expect hydrody- 
namics to play an important role in their dynamics. In 
particular, as the defects move, the coupling between the 
changing director field and the velocity field (so-called 
back-flow) may play a signiflcant part in the motion. Ex- 
perimental evidence [Q shows that this is indeed the case 
when a nucleated domain where the director field is hor- 
izontal grows in a twist or vertical environment due to 
the influence of the surfaces or an external electric field. 
The speed of the domain boundary is found to depend 
strongly on the local defect configuration. 

Previous investigations of defect dynamics have 

either ignored the flow field or taken account of its effect 
phenomenologically. Here we aim to generalize this work 
by treating the full hydrodynamic equations of motion 
for a nematic liquid crystal. We consider the annihila- 
tion of a pair of defects of strength s = ±1/2. We find 
that back-flow can change the speed of defects by up to 
~ 100%. Defects of different strength couple to the flow 
field in different ways. This leads to a dependence of 
speed on strength which can occur either as a result of 
the back-flow or if the elastic energy is treated beyond a 
one-elastic constant approximation. 

The hydrodynamics of liquid crystals is often well de- 
scribed by the Eriksen-Leslie-Parodi equations of motion, 
which are written in terms of the director field. How- 
ever these are restricted to an uniaxial order parameter 



of constant magnitude. Thus they are inadequate to ex- 
plore the hydrodynamics of topological defects where the 
magnitude of the order parameter has a steep gradient 
and becomes biaxial within the core region Here we 
consider the more general Beris-Edwards [|o| formulation 
of nematohydrodynamics, where the equations of motion 
are written in terms of a tensor order parameter Q. 

The equilibrium properties of the liquid crystal are de- 
scribed by a Landau-de Gennes free energy density [1). 
This comprises a bulk term 
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which describes a first order transition from the isotropic 
to the nematic phase at 7 = 2.7, together with an elastic 
contribution 
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where the L's are material specific elastic constants. The 
Frank expression for the elastic energy, written in terms 
of the derivatives of the director, can be simply mapped 
to (H) |]lo|. A controls the relative magnitude of ft and 
fd- The Greek indices label the Cartesian components of 
Q, with the usual sum over repeated indices. 

The dynamics of the order parameter is described by 
the equation 
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where u is the bulk fluid velocity and F is a collective ro- 
tational diffusion constant. The term on the right-hand 
side of equation (||) describes the relaxation of the order 
parameter towards the minimum of the free energy J- 
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The term on the left-hand side is 



S(W, Q) = (^D + n)(Q + 1/3) + (Q + I/3)(eD - U) 
-2e(Q + I/3)Tr(QW), (5) 
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where D = (W + W^) /2 and n = (W ~ W^) /2 are the 
symnietric part and the anti-symmetric part respectively 
of the velocity gradient tensor Wa/3 = dpUa- £, is related 
to the aspect ratio of the molecules. 

The velocity field u obeys the continuity equation and 
a Navier-Stokes equation with a stress tensor generalized 
to describe the flow of nematic liquid crystals 
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where p and T are the density and temperature. Notice 
that the stress (^) depends on the molecular field H and 
on Q. This is the origin of back- flow. Details of the 
equations of motion can be found in reference |ic| ]. The 
equations (|l|)-@ were solved numerically using a lattice 
Boltzmann algorithm described in [pT[ . 

Consider a pair of defects of topological strength s = 
±1/2 situated a distance D apart in a nematic liquid 
crystal, as shown in Fig. [|(a). We consider a two- 
dimensional cross section of the two line defects, assum- 
ing that the order parameter does not change in the per- 
pendicular direction (although the director may point out 
of this plane) . The two defects are topologically distinct 
only in two dimensions, but even in three dimensions they 
are separated by an energy barrier for the typical elastic 
constants we study here. 

A phenomenological equation of motion can be written 
down by assuming that the attractive force between the 
two defects [|j is counterbalanced by a friction force Q 
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where D is the defect separation, Rc is the defect core 
size and po is a constant. 

In Ref. the director field and the trajectory of the 
defects were obtained analytically and the defect veloc- 
ity was determined as a function of parameters of the 
medium. (A review of the earlier development of the 
theory of defect dynamics is also given in [^.) However 
this and, as far as we are aware, all other analyses of de- 
fect dynamics have ignored back-flow. This means that 
the approach to equilibrium is relaxational, determined 
entirely by the derivative of the free energy with respect 
to the order parameter, with the flow playing no role. A 
further simplification in previous work is that the Frank 
elastic constants were assumed to be equal. 

We can examine relaxational dynamics using a 
Ginzburg-Landau equation for the director field [p^, i.e. 
Eqn. with the velocity set to zero. The Ginzburg- 
Landau equation with a single elastic constant is invari- 
ant under a local coordinate transformation mirroring 



the director on the x axis (where we define x as the axis 
connecting the two defects cores). This corresponds to 
the transformation 
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The order parameter fields of the two defects with topo- 
logical charges s = ±1/2 transform into each other. Thus 
approaches based on a simple Ginzburg-Landau equation 
predict that when the defects move they follow symmet- 
ric dynamical trajectories. 

Fig. ^a) shows the position of two annihilating topo- 
logical defects, with topological charge s — +1/2 (upper 
curve) and s = —1/2 (lower curve), as a function of time. 
Fig. ^b) shows the velocities of the defects as the func- 
tion of their separation D. (See footnote for the 
simulation parameters.) 
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FIG. 1. (a) The director field of two annihilating topolog- 
ical defects with strength s = ±1/2. (b) Velocity field of the 
two defects. 

Consider first the dashed trajectories. These were 
obtained with the flow field switched off, the case for 
which equation (^ reduces to a Ginzburg-Landau model. 
As expected the two defects move with the same speed 
and annihilate halfway between their initial positions 
(marked by a horizontal line in Fig. ^(a)). The results for 
the Ginzburg-Landau model fit the simple formula (^ for 
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/io = 124/im2/s and Rc = 0.0233^m for D > R^. Around 
D ^ Rc the formula overestimates the defect speed. 
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FIG. 2. (a) Positions of the defects as a function of time. 
Upper (lower) curves are for s — +l/2( — 1/2). The different 
lines correspond to a Ginzburg-Landau model(dashed), the 
hydrodynamic equations of motion with F — 6.25 Pa~^s~^ 
(solid) and with F = 7.76 Pa"^s~^(dotted). (b) Defect 
speed as a function of separation. Upper and lower solid 
curves: the s = +1/2 and s — —1/2 defect trajecto- 
ries with hydrodynamics and F = 6.25; Dashed curve: 
Ginzburg-Landau model. The -1-1/2 defect is considerably 
(100%) accelerated for D > 0.25/im compared to the results 
of the Ginzburg-Landau model. The speed of the —1/2 defect 
is only slightly affected by the back-flow. 

Consider the effect of the back-flow under the trans- 
formation of Eq. . Examining the stress tensor one can 
see that the last term in Eq.(^) does not change under 
the transformation but the off-diagonal elements of the 
other terms have their sign inverted. This reflects the 
two sources of back-flow in this problem. The first has 
to do with the defect core. Order is suppressed in the 
core, which results in an increase in viscosity at the core 
(The isotropic viscosity a4, is proportional to (1 — g)^ 
where q is the magnitude of the order [^^). As a re- 
sult, the movement of the core induces similar vortices to 
those produced when a solid cylinder is moved through 
a fluid. This flow is independent of the sign of the de- 
fect, and points into the direction of defect propagation 
at the core. The second source of back-flow comes from 
the reorientation of the director field away from the core. 



This flow depends on the sign of the spatial derivative of 
the director orientation. As a result, these two sources 
of flow reinforce in one case and partially cancel in the 
other, thus giving the anisotropy. The flow fleld around 
the defects is depicted in Fig. |l|(b). A strong velocity 
vortex pair is formed around the +1/2 defect, with the 
flow pointing in the direction of defect propagation. The 
flow around the —1/2 defect is much weaker, and points 
opposite to the direction of defect motion. 

The solid line in Fig. |^(a) corresponds to a simulation 
of the full hydrodynamic equations of motion 
There is a marked decrease in the time to coalescence. 
This is primarily because the speed of the s = +1/2 de- 
fect is increased by ^ 100% compared to the case with- 
out flow for D > 0.25/im. The s = — 1/2 defect is only 
slightly affected by the back-flow; the change in veloc- 
ity is less than 20%. Due to the speed anisotropy the 
defects do not meet halfway between their initial posi- 
tions. For defects initially 1/im apart the displacement 
of the coalescence point (Axi = 0.149/im) is smaller 
than might be expected from the substantial speed-up 
of the s = +1/2 defect. This is because the relative flow- 
induced increase in velocity drops dramatically near the 
defect core {D < 0.25/xm) where the defects are moving 
the fastest. At these short separations the relaxational 
dynamics dominates the hydrodynamics. 

Changing the various material parameters of the sam- 
ple will affect the velocity of the defects and their speed 
anisotropy. We find, as expected, that as the viscosity 
in the Navier-Stokes equation increases the motion ap- 
proaches that of the Ginzburg-Landau model. Increas- 
ing r in increases the speed of relaxation to the min- 
imum of the free energy, thus increasing the speed of 
the defects. The speed anisotropy however decreases be- 
cause the weight of the free energy relaxation process is 
increased relative to the hydrodynamics. For example 
the dotted curve in Fig. ^(a) corresponds to F = 7.76 
Pa"'^ s~^ (compared to F = 6.25 Pa~^ s"'^ for the solid 
curve). The displacement of the coalescence point is 
Ax2 = 0.128/im < Axi. Decreasing A in the free energy 
(|l|) increases the defect size. The defects move faster but 
the velocity disparity decreases, again because the impor- 
tance of the relaxational dynamics is increased relative to 
the hydrodynamics. 

The results in Fig. || are for a single elastic constant 
(equal Frank elastic constants or, equivalently, L2 = 
L3 = 0). If a more general model for the elasticity is 
considered, allowing L2 or L3 ^ 0, the invariance of 
the Ginzburg-Landau equation under the transformation 
(^) is broken. Therefore one might expect a difference in 
the flow velocities of s = ±1/2 defects even if back-flow 
is not considered. 

This is indeed the case. If L2 = and L3 < (L3 > 0) 
the s — +1/2 (s = —1/2) defect moves faster. For ex- 
ample, if Li = 8.73 pN, L2 = 0, and L3 = 15.88 pN 
Q a comparison of the velocities f -1-1/2, t^-i/2 of the 
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s = 1/2 and s = — 1/2 defects respectively gives a speed 

anisotropy a„ = (u+1/2 - t^-i/2)/( "^'''%""'^" ) ~ -13% 
at a defect separation D = 0.5/im. For comparison, the 
anisotropy caused by the back-flow is ~ +68%. The dis- 
placement of the coalescence point is Aa: = 0.029/im. For 
most materials L3 > (corresponding to Kn < K33), 
leading to a speed anisotropy opposite to that arising 
from the hydrodynamics. 

If L2 7^ and L3 = the velocity anisotropy is very 
small since for these values of the elastic constants the 
relaxational Ginzburg-Landau dynamics remains invari- 
ant under the mapping in the limit of uniaxiality and 
constant magnitude of the order parameter [ p^ . The 
speed anisotropy is small because these conditions are 
only relaxed within a defect core. For example, L3 — 
and L2 = 15.88 pN @ leads to |a„| < 2%. 

Anisotropics in the speed of domain walls of up to 50% 
have been observed in experiments on pi-cell liquid crys- 
tal devices where the movement of twist and splay-bend 
walls is important in mediating the formation of the oper- 
ating (bend) state from the ground (splay) state. Defects 
form spontaneously at these walls and preliminary simu- 
lations show that back-flow effects are responsible for the 
velocity anisotropy |T^. It is also of interest to investi- 
gate the role of defect motion in many other new gener- 
ation liquid crystal devices. For example multi-domain 
nematic modes improve viewing angles at the expense of 
introducing defects into the director profile and under- 
standing the behavior of such defects as the electric field 
is varied will help control device performance. In zenithal 
bistable nematic devices switching is between two (meta) 
stable zero-field states. Switching between the states is 
mediated by the movement of topological defects (is) . 

To conclude, we have used a formulation of nemato- 
dynamics based on the tensor order parameter to study 
the hydrodynamics of topological defects in nematic liq- 
uid crystals. We find that the coupling between the order 
parameter field and the flow has a significant effect on de- 
fect motion: in particular it introduces a substantial dif- 
ference between the velocities of defects of different topo- 
logical charge. Similar but smaller velocity anisotropics 
can result from changing the elastic constants. 
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